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RATES OF OBSERVABLE BLACK HOLE EMERGENCE IN SUPERNOVAE 
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ABSTRACT 

A newly formed black hole may be directly identified if late-time accretion of material from the base of the 
ejected envelope generates a luminosity that is observable in the tail of the supernova light curve. In this work we 
estimate the rate at which events where the black hole "emerges" in the supernova light curve can be detected with 
present capabilities. Our investigation is based on an analytical model of the accretion luminosity at emergence 
as a function of progenitor mass, coupled to the inferred rate of observed Type II supernovae in nearby galaxies. 
We find through a parameter survey that under optimistic assumptions the potential rate of observable events can 
be as high as several per year. However, supernovae which produce black holes are also likely to be low energy 
explosions and therefore subluminous, as was the case for the best candidate to date, SN1997D. If black hole- 
forming supernovae are underdetected owing to lower luminosities, the rate of observing black hole emergence is 
probably not larger than once every few years. We therefore emphasize the importance of dedicated searches for 
nearby supernovae as well as faint supernovae projects for improving the prospects of observationally certifying 
the supernova-black hole connection. 

Subject headings: accretion, accretion disks — black holes — supernovae: general 



1. INTRODUCTION 

Theoretical studies of core-collapse supernovae suggest that 
the remnant they leave behind can be either a neutron star or a 
black hole. The nature of the remnant is determined mostly by 
the amount of material the collapsed core accretes while the ex- 
plosion is still progressing: if the final mass of the core exceeds 
the maximum mass it can sustain, it will collapse to a black hole 
rather than evolve to a stable neutron star. 

Observational support for this theoretical scenario is not 
complete. Radio pulsar emission and x-ray observations pro- 
vide observational evidence associating neutron stars with sites 
of known supernovae, but similar evidence for a black hole - 
supernova connection is still mostly unavailable. Indirect ev- 
idence that the black hole candidate in the x-ray binary sys- 
tem GRO J1655-40 was formed in a supernova explosion has 
recently been reported by Israelian et al. (1999). They base 
their inference on the relatively high abundances of nitrogen 
and oxygen on the surface of the companion, which has too 
low a mass to have produced such abundances by thermonu- 
clear burning. Assuming those elements were deposited on the 
companion by the supernova explosion of the primary star pro- 
vides indirect evidence of a black hole formed in a supernova. 

In a series of recent papers we examined the prospects of di- 
rectly identifying a black hole created in a supernova through 
late-time accretion onto it (Zampieri et al. 1998a,b; Balberg et 
al. 2000). Such accretion is powered by material from the base 
of the ejected envelope which remains bound to the black hole 
and gradually falls back and is accreted. Analytic (Colpi et 
al. 1996) and numerical (Zampieri et al. 1998a) investigations 
showed that the rate of late-time spherical accretion, when it 
becomes dust like, is expected to decline as a power law in time 
and so is the corresponding accretion luminosity. Such an ac- 
cretion luminosity will produce a distinct signature on the total 



light curve if and when it becomes comparable to the output 
of other power sources, i.e., the initial internal energy of the 
envelope and decays of radioactive isotopes synthesized in the 
explosion. If detected, the unique time dependence of the ac- 
cretion luminosity due to fall back would be a direct sign of the 
black hole "emerging" in the supernova light curve. 

To date, black hole emergence in supernovae is yet to be ob- 
served. The prominant obstacle is that radioactive heating in the 
light curves of typical core collapse supernovae obscures any 
accretion luminosity for very long times. Even though radioac- 
tive heating decays exponentially while accretion luminosity 
declines only as a power law, the absolute power in radioac- 
tive heating remains the dominant source of the light curve un- 
til the supernova is no longer detectable. The case of SN1987A 
is a useful example (Zampieri et al. 1998a): while no firm ev- 
idence of a newly formed neutron star has been found, posi- 
tive indication of the presence of a black hole (if one was pro- 
duced) would be unobservable due to radioactive heating for 
about 1000 years. 

Viable candidates for observing black hole emergence are su- 
pernovae which show very little radioactive heating. In princi- 
ple, larger mass remnants, which are more likely to be black 
holes (Fryer 1999), are also more likely to be associated with 
supernovae which show diminished abundances of radioactive 
isotopes. Radioactive isotopes are produced in the deepest 
part of the exploding star, and therefore can be absorbed in a 
larger remnant in the earliest stages of fallback (Woosley & 
Weaver 1995). The best two recent candidates for supernovae 
which ejected very little radioactive isotopes and possibly pro- 
duced black holes were SN1994W and SN1997D. In SN1997D 
the abundance of 56 Co inferred in the light curve was only 
2 x 10 -3 Mq, and the explosion energy was estimated as only 
4 x 10 50 ergs, both significantly lower than typical in Type II su- 



pernovae (Turatto et al. 1998). Our previous analyses (Zampieri 
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et al. 1998b; Balberg et al. 2000) focused on SN1997D and 
showed that a 3 Mq black hole in its remnant may emerge 
about 1000-1200 days after the explosion. Unfortunately, even 
at emergence the luminosity would have been only marginally 
detectable with HST. In the case of SN1994W 2 x 10~ 3 M of 
56 Co was the upper limit for radioactive heating, so the accre- 
tion luminosity might have emerged much earlier (and at higher 
value than in SN1997D). Unfortunately (again) the tail of the 
light curve was apparently dominated by circumstellar interac- 
tion (which was absent in SN1997D), and thus obscured any 
weaker source of heating from inside the envelope (Sollerman 
etal. 1998). 

Can we expect better success in a systematic quest for ob- 
serving black hole emergence in supernovae? Our goal here is 
to perform a first estimate of the potential rate of observations, 
and examine possible strategies for improving the prospects of 
such an observational search. We base our estimate on an ana- 
lytic approximation for the luminosity at emergence of a black 
hole in a supernova as a function of progenitor mass, and then 
couple it to estimates of the rate of supernovae in nearby galax- 
ies. Although the model is crude and excludes several details 
which determine specifics of supernovae explosions, their light 
curves and remnant formation, we are able to point to some un- 
derlying features. Our conclusions suggest that the likely rate at 
which black hole emergence can be observed with present ca- 
pabilities and supernovae searches is probably once every sev- 
eral years, although with favorable assumptions this rate can be 
as high as about once a year. We also find that the dominant 
factor in determining this rate is most likely the difficulty of de- 
tecting nearby faint supernovae. A systematic search for such 
supernovae would significantly enhance the probability of dis- 
covering black hole emergence and observationally confirming 
the black hole - supernova connection. 

We begin in § 2 where we develop the model that relates the 
emergence luminosity to progenitor mass. This function serves 
as the basis for estimating the values of emergence luminosi- 
ties, presented in § 3. The main results are given in § 4, where 
the emergence luminosities are coupled to a distribution of pro- 
genitor masses, allowing for estimates of the rate of observable 
supernovae in nearby galaxies and providing us with an esti- 
mate of the rate of observable black hole emergence events. 
Conclusions and discussion are given in § 5. 



2. THE MODEL 

To assess the observability of black hole emergence we must 
first determine the accretion luminosity of the black hole and its 
value at emergence. We employ a simplified model in which the 
emergence luminosity, Lbh, is a unique function of progenitor 
mass, . This is clearly a crude approximation, since explo- 
sion energy, black hole mass and post-shock evolution are likely 
to be dependent on other factors, such as progenitor rotation 
(Fryer & Heger 2000), metallicity (Woosley & Weaver 1995) 
and effects of a binary companion (Willstein & Langer 1999). 
Nonetheless, such an approach is a useful first approximation 
for deriving the observable event rate (and indeed is frequently 
used in other applications regarding core collapse supernovae, 
such as determining the neutron star and black hole initial mass 
functions; see Timmes, Woosley & Weaver 1996, Fryer 1999, 
Fryer & Kalogera 1999). We discuss some of the limitations 
and uncertainties concerning the unique function Lbh (M* ) to- 
wards the end of this section. 



2.1. The Accretion Luminosity 

Constructing the function Lbh(M*) requires the time- 
dependence of the luminosity due to accretion onto the black 
hole. As shown analytically by Colpi et al. (1996) and numer- 
ically by Zampieri et al. (1998a), at late times after the explo- 
sion, spherical accretion onto the black hole becomes dust-like, 
and the accretion rate, M, declines asMa i~ 5//3 . Zampieri et 
al. (1998a) also showed that the resulting accretion luminosity 
is consistent with the Blondin (1986) estimate of hypercritical 
accretion, where L acc oc M 5 ' 6 , so that 

L acc cx t- 25 / 18 . (1) 
This is a fundamental prediction of spherical accretion onto a 
black hole following a supernova, and can serve as distinct sig- 
nature of accretion in the overall light curve (Zampieri et al. 
1998a; Balberg et al. 2000). We address the question of likeli- 
hood of spherical accretion later in this section. 

The proportionality constant in equation (1) is dependent on 
the black hole mass, M bh, and on the properties of the bound 
material, including composition, density and velocity profiles. 
For the purpose of an analytic estimate an effective time, tdust, 
can be defined so that 
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where LEdd = AkcGMbh / 'k is the Eddington luminosity for 
material with opacity k accreting on a black hole of mass Mbh- 
For bound material which initially has a uniform density, p n , 
and a homologous velocity profile, v(r, t = 0) = r/t n (where 
to is the expansion time scale), the time t^ust can be roughly 
estimated as (Balberg et al. 2000, eq. [30]) 
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10 9 gm cm~ 3 s 

where fj, is the mean molecular weight of the bound material. 
This estimate for tdust includes the effect of radiation pres- 
sure during the earliest phase of fallback, when an accretion 
luminosity at about the Eddington limit modifies the flow of 
the bound material. This seems to be the case for realistic su- 
pernova envelopes (see Balberg et al. 2000 for details). We 
emphasize that the time tdust is an effective quantity: at that 
actual time after the explosion the accretion is not expected to 
have settled into dust-like motion. Rather, tdust corresponds to 
the time when the luminosity extrapolated back from late-time 
equals the Eddington limit. 

2.1.1. Explosion Energy and Black Hole Mass 

Equation (3) introduces the quantities required for estimating 
the accretion luminosity. For a given progenitor, these quanti- 
ties are determined primarily by the explosion energy and the 
manner in which it is distributed in the envelope. We follow 
the approach of Fryer & Kalogera (1999), and parameterize the 
explosion by the total energy available for the through core col- 
lapse, E exp , and the fraction, /, of this energy which is spent 
on unbinding the star. The remainder of the explosion energy is 
converted mostly to kinetic energy of the ejected material. 

The energy available to unbind the star can be combined with 
the binding energy profile of the progenitor to calculate the rem- 
nant mass, M rem , as a function of M* according to 
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FIG. 1. — Estimate of the total explosion energy, E exp (eq. [5]) and pre- 
dicted remnant mass, M rem as a function of progenitor mass, M«, for dif- 
ferent values of the parameter / (top to bottom: / = 0.25, 0.50 and 0.75). 
Data points correspond to the values from the numerical model (eq. [4), lines 
are linear best fits for the M rem (M*) function for high mass remnants. The 
solid line at M rem = 3 Mq marks the assumed threshold mass for black 
hole formation. 

where EBE{m) is the binding energy profile, for which we use 
the massive star models of Woosley & Weaver (1995). Figure 1 
shows the remnant mass versus progenitor mass relation found 
with / = 0.25, 0.5 and 0.75 where we adopted an analytic ex- 
plosion energy function of 
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ergs , (5) 



for M* in units of Mq. This energy function (also plotted in 
Fig. 1) is based on a fit to 2D core-collapse simulations of pro- 
genitors in the range 8 < M* < 40 Mq, evaluated as the change 
in energy before collapse and after about 1 second of the ma- 
terial that experiences explosive shock heating. (Fryer 1999). 
Note that for the higher range of progenitor masses, which are 
the candidates for forming black holes, the remnant mass varies 
about linearly with the progenitor mass. Obviously, the remnant 
mass increases with decreasing /. 

The nature of the remnant depends on whether its mass ex- 
ceeds the maximum mass of a neutron star. As a conservative 
lower limit we will require M rem > 3 M Q (baryonic mass, 
which is that found in eq. [4]) for the remnant to be a black 
hole. This consistent is with most realistic equations of state 
(see Cook, Shapiro & Teukolsky 1994). 

2. 1 .2. Properties of the Bound Material 

The estimate of the accretion luminosity includes a depen- 
dence on both the composition and kinetic energy of the bound 
material. For a homologously expanding medium the combi- 
nation po^o can b e expressed in terms of the global quantities 
of the post-shock helium-rich layer, which is the source of the 
late-time accretion (Balberg et al. 2000): 

p (Re)4 - 0.222 x 10 9 (^) 5/2 (6) 

-3/2 



*WHe) \ ' 3 

io«erg ) gm cm s 



where M# e is the mass of the layer and E^ in (He) is its total 
kinetic energy. 

The mass of the helium layer can be estimated by subtract- 
ing the remnant mass calculated in equation (4) from the mass 
of the helium core of the progenitor. We parameterize the ki- 
netic energy of the helium layer as a fraction, /h , of the total 
kinetic energy delivered to the ejected envelope: E kin (Yic) = 
/He(l — f)E exp . Simulations of Type II explosions typically 
suggest that /h c is of the order of a few percent (Nomoto et al. 
1994; Arnett 1996), generally decreasing for larger mass pro- 
genitors. The exact value does, however, depend on the details 
of the explosion and on progenitor structure. 

2.2. Emergence Time and Luminosity 

The black hole mass and the effective time tdust uniquely de- 
termine the time dependent accretion luminosity (eq. [2]) which 
can then be compared to other sources of emission in the super- 
nova. Specifically, we define the time of emergence, tsH, as 
that when the accretion luminosity reaches | of the total bolo- 
metric luminosity. This time and the corresponding emergence 
luminosity, Lbh = L acc (tBH), depend on the source of lumi- 
nosity with which accretion competes. This source is either the 
initial internal energy deposited in the envelope by the explo- 
sion, or ongoing heating due to decay of radioactive isotopes. 

The radioactive isotopes relevant to powering the light curve, 
56 Ni and its decay product 56 Co , 57 Co and 44 Ti , are ex- 
pected to form in the deepest layers of the supernova enve- 
lope (Timmes et al. 1996). Hence, their final abundances in 
the ejected envelope may be depleted with respect to the to- 
tal quantities produced if a significant fraction of these inner 
layers falls back onto the collapsed core while the explosion 
is still progressing. In the simplest picture of black hole for- 
mation in an otherwise "successful" supernova, larger remnant 
masses coincide with more fall back and therefore with lower 
abundances of ejected radioactive isotopes. A quantitative man- 
ifestation of this assessment is found in survey by Woosley & 
Weaver (1995). For a wide range of parameters, they find that 
explosions which produce remnants with M rem > 3 M Q eject 
negligible abundances of 56 Ni , while in explosions which leave 
behind remnants of 1 — 2.5 Mq, the amount of ejected 56 Ni is 
of order 0.1 M Q (the exact value depends on the details of the 
progenitor and explosion). 

2.2.1. No Radioactive Isotopes 

If we adopt this simple picture of radioactive isotope pro- 
duction in supernovae, we may assume that supernovae which 
produce black holes eject negligible amounts of radioactive iso- 
topes. Correspondingly, the accretion luminosity must compete 
only with the thermal emission of the internal energy deposited 
during the explosion. 

Immediately after the explosion the envelope is highly ion- 
ized, and therefore opaque to its own thermal photons. The 
internal energy is emitted slowly through photon diffusion, and 
the expansion of the envelope is roughly adiabatic, with average 
density and average temperature falling off inversely with time 
(Arnett 1996; Zampieri et al. 1998a). Eventually, expansion de- 
grades the envelope temperature sufficiently to allow its mate- 
rial to recombine. The recombined material is practically trans- 
parent to the thermal photons, so the envelope cools rapidly. 
The process is so rapid that the photon distribution cannot ad- 
just to recombination - rather, a recombination "front" sweeps 
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through the envelope (Arnett 1996), liberating most of its ther- 
mal energy. 

Recombination sets in after a recombination time 



tree — io(H) 



Jo(H) 
T rec (H) 



(7) 



where T (H) is the initial average temperature of the hydrogen 
envelope (which contains most of the mass and the internal en- 
ergy), io(H) is its initial expansion time and T rec (H) is a repre- 
sentative temperature at which hydrogen recombines. Roughly, 
the internal energy is emitted over a time ~ t rec , during which 
the recombination front sweeps recedes through the bulk of the 
envelope. We therefore assume that the accretion luminosity 
emerges at a time of roughly 



(8) 



This estimate is consistent with observed Type II supernovae 
(SNell) light curves, where the time until the recombination 
peak is roughly equal to the time it lasts until dropping below 
luminosity from 56 Co emission. For example, in SN1987A the 
recombination peak dominated the light curve between about 
50 and 120 days after the explosion; in SN1997D, it appears 
that the recombination peak lasted between about 40 and 80 
days after the explosion. This estimate is also consistent with 
numerical simulations of SNell light curves, and, specifically, 
with our previous simulations of black hole emergence in the 
absence of radioactive heating (Zampieri et al. 1998a; Balberg 
etal. 2000). 

In the context of our analytic model, we may estimate the 
recombination time by approximating the entire ejected ma- 
terial (both the helium-rich layer and the hydrogen-rich enve- 
lope) as uniform both in temperature and in expansion time, 
so To(H) = To, to(H) = to. Both quantities can then be 
expressed in terms of the initial properties of the progenitor, 
namely the mass of the hydrogen envelope, Mh, and the initial 
outer radius, R : 
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(9) 



(10) 



In equations (9-10), E t h is the total initial thermal energy in 
the envelope, while Vo(H) and £^„(H) are the outer veloc- 
ity and the total kinetic energy of the hydrogen-rich envelope, 
respectively. Now Ekin (H) is defined by our previous parame- 
terization, 



Ekin (H) = Ekin — Ekm (He) — 

(1 - fuc)Ekin = (1 — /Hc)(l - f)E exp 



(11) 



The thermal energy, E t h, is some fraction of the kinetic en- 
ergy, E t h = fthEkin- To be consistent with the assumption 
of homologous ballistic expansion, the thermal energy must be 
significantly lower than the kinetic energy, i.e., f t h < 1; in 
the following we will use the value f t h = 0.1. Note that this 
choice does not necessarily reflect the actual conditions in the 
supernova at shock break out, when the external radius is Rq, 
since the envelope will not yet have settled into a homologous 
profile. Rather, it is an effective parameterization of the post- 
shock flow once it settles into homologous expansion (after a 
few expansion times; Arnett 1996). 



2.2.2. Some Radioactive Isotopes 

Even relatively low abundances of 56 Co , 57 Co and 44 Ti will 
lead to the tail of the light curve being dominated by radioac- 
tive heating. In the case of SN1997D the observed amount of 
2 x 10~ 3 Mq was found to be sufficient to swamp the accretion 
luminosity for about 1000 days after the explosion (Zampieri et 
al. 1998b; Balberg et al. 2000). 

It is difficult to generalize the significance of the observed 
56 Co abundance in SN1997D, since there have been only a 
handful of other supernovae where the 56 Co ( 56 Ni ) abundance 
was observed to be below 0.01 Mq (Young et al. 1998), with 
SN1997D being the only case where this abundance has ac- 
tually been determined. In view of the importance of radioac- 
tive heating to the prospects of observing black hole emergence, 
we consider several models of the relation between progenitor 
mass and abundance of radioactive isotopes. We quantify these 
models through the abundance of 56 Co , which dominates ra- 
dioactive heating after recombination ( 57 Co and 44 Ti are dis- 
cussed below). The black hole emerges in the light curve 
when the accretion luminosity equals the luminosity due to 
56 Co heating, L acc (t BH ) = Lse Co (t BH ). 

As an extreme model, we examine the possibility that due to 
mixing in the earliest stages of the explosion a minute amount 
of radioactive isotopes always reaches the layers which are 
eventually ejected, regardless of the final black hole mass. 
Hence we assume that all black hole forming supernovae eject 
M56 Co = 2 x 10~ 3 M Q , as was inferred for SN1997D. As 
a compromising alternative, we also consider a model where 
there is always finite 56 Co abundance which does depend on 
the black hole mass: for lack of an actual physical model, we 
stipulate that M5e Co = 2 x 10~ 3 M© for M B h = 3 Mq, but 
declines rapidly with increasing black hole mass through the 
relation 



M 5 6 Co {M BH ) =M exp 



M, 



BH 



(12) 



,3 M 0/ 

assuming M BB > 3 Mq. We consider a "moderate de- 
cline" of (3 = 2 and a "rapid-decline" of /? = 8. We chose 
M = 5.44 x 10~ 3 Mq, so that M™ Co (M BH = 3 Mq) = 
2 x 10~ 3 Mq. 

At times greater than a few tens of days (when 56 Ni decay is 
complete) the luminosity due to heating from 56 Co decays can 
be estimated as (Woosley et al. 1989) 

'M5e Co " 



L56Co{t) 
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r56 f 



(13) 
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where e 7 ( 56 Co) 

-|41 



[e 7 ( 56 Co).f 7 (t)+e e+ ( 56 Co)]e 

1.27 x 10 43 ergs s^and e e+ ( 56 Co) = 
4.45 x 10 4i ergs s _1 are the energy emission rate in 7— rays 
and positrons (produced per 1 Mq of 56 Co ), respectively, and 
t( 56 Co) = 111.3 days is the life time of 56 Co . The coefficient 
/ 7 includes the effects of a finite optical depth of the envelope 
to the emitted 7— energy photons. Photons which escape before 
thermalizing their energy do not contribute to the bolometric 
light curve, i.e., / 7 < 1. Roughly, / 7 can be estimated through 
the time dependent optical depth of a layer with width R(t): 



1 - fiit) w exp[ 
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(14) 



where the opacity of the material to the 56 Co decay pho- 
tons; a good estimate is that for the typical energy of these pho- 
tons, k 7 = 0.03(1 + X) cm 2 gm _1 , where X is the hydrogen 
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mass fraction of the material. For a uniform density layer with a 
mass M' and kinetic energy E' kin the 7— energy photon optical 
depth can be estimated as (Balberg et al. 2000): 



R(t) 

I 



K 7 p(r)dr w K~fR out (t)p(t) = k 



3 3 M 



(15) 



and we used R out (t) = V a t (V a is the velocity at the outer 
edge of the layer) and equation (6). Note that the relation 
/ 7 cx 1 — cxp(— 1 / 1 2 ) combines with the natural exponential to 
accelerate the temporal decline of Lseco (Eq. [13]), much faster 
than the power-law decline of L acc (Eq. [2]). 

2.3. Summary: Model Parameters and Uncertainties 

The parameterizations introduced above allow to determine 
the luminosity at emergence as a function of progenitor mass. 
The key quantity is the effective time tdust '■ 
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In the following we assume that the composition of the bound 
material as a mixture in mass of 0. 1 hydrogen, 0.45 helium and 
0.45 oxygen, which is typical of SNell models (T. R. Young, 
1999, private communication). For this composition [i = 1.26 
and k = 0.22. 

In the case of no radioactive isotopes in the envelope, the 
emergence time is determined by recombination, and we have 
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In the case of a finite 56 Co abundance, the luminosity due 
to radioactive heating is estimated with equations (12-13). 
The 7— ray optical depth due to the helium-rich layer and the 
hydrogen-rich envelope is evaluated using equations (14-15). 
For SNell we assume that in the helium-rich layer X = 0.1 
and that in hydrogen-rich envelope X = 0.7, giving rise to the 
following result for the 7— ray trapping factor: 
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For most of realistic values of SNell, the contribution of the 
helium-rich layer (first term in eq. [18]) dominates 7— ray trap- 
ping, since Mhc/Mq > Mh /IQMq. However, in the low 
energy explosions of the most massive progenitors the contri- 
bution of the hydrogen-rich envelope becomes comparable, and 
must be included. 



2.4. Comments on the Model 

The analytic form of equations (16-18) enables us to assess 
the time of emergence and the corresponding luminosity for a 
well defined set of assumptions. Naturally, this convenience is 
achieved at the price of several significant simplifications, and 
we note some obvious sources of uncertainties. 

Homologous Expansion. The assumption of a homologous 
expansion profile is in itself a crude approximation (see the 
comparison of an analytic estimate for the emergence time and 
luminosity with the results of a simulation with a nonhomol- 
ogous initial profile in Balberg et al. 2000). Furthermore, su- 
pernova simulations (see, e.g., Woosley 1988, Shigiyama & 
Nomoto 1990) also suggest that the helium-rich layer and the 
hydrogen-rich envelope do not share a common expansion time, 
but rather t (H) < t (He). The ratio of the two tends to be 
very model dependent. Our model must therefore be treated as 
illustrative, rather than exact. 

Secondary Model Simplifications. The lower mass limit for 
black holes may be significantly less than 3 M Q (Bethe & 
Brown 1995). In principle this could add many more events, 
but since lower mass remnants (regardless of nature) are likely 
to be correlated with larger amounts of radioactive isotopes in 
the ejecta, they would not make favorable candidates for ob- 
serving black emergence. Also, there is a small inconsistency 
since the bound material for late-time accretion is actually in- 
cluded in Mbh in the recipe of equation (4). However, since 
Mbh > 3 Mq while the amount of bound material for late- 
time fallback is 0.1 — 0.2 Mq, the impact of this inconsistency 
is small. Finally, we note that the estimate of t rec in equation 
(17) includes two additional parameters, namely the effective 
thermal energy fraction and the progenitor radius, making the 
derivation less general. However, since the dependence on both 
is weak, we do not introduce significant errors by using only 
representative values for these parameters. 

Explosion Energy. Our one parameter model of Lbh(M^) 
is obviously a simplification, since in reality there should be 
a range of explosion energies for any progenitor mass. Such 
a range would be due to rotation, metallicity and other fac- 
tors which affect progenitor structure and evolution. We aim 
to accommodate this simplification by conducting a parameter 
survey with the values of / and fn e , which should span a rea- 
sonable range of uncertainty in this regard. 

Mass Loss. The model does not include potential mass loss, 
which is expected for massive stars. Mass loss may be impor- 
tant in determining the eventual light curve of the supernova, 
and perhaps even the properties of the remnant. In this context 
we distinguish between modest mass loss expected in single 
stars and the substantial mass loss that occurs in a binary sys- 
tem. In the former the star is likely to retain a sizable fraction 
of its hydrogen envelope, and so explode as SNII (Willstein & 
Langer 1999) without a significant effect on the explosion en- 
ergy and the black hole mass (MacFadyen et al. 1999; Fryer 
& Kalogera 1999). The more efficient mass loss from massive 
stars in binaries is actually likely to strip the hydrogen envelope 
from the progenitors (Willstein & Langer 1999), and if it occurs 
prior to helium ignition, it may eventually even affect the nature 
of the explosion and remnant (Fryer & Kalogera 1999), compli- 
cating a simple estimate along the lines of our model. However, 
once the hydrogen envelope is stripped from the star, the likely 
result will be a Type Ib/Ic supernova, which is unfavorable for 
detecting black hole emergence (see below). Hence, we can 
avoid considering the effect of significant mass loss on explo- 
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sion properties and results. On the other hand, modest mass 
loss does increase the likelihood that circumstellar interaction 
and dust formation will obscure the emergence luminosity, and 
therefore may have an unfavorable effect on the prospect of de- 
tecting emergence. We return to this point later in § 4. 

Spherical Accretion. We must caution that the results re- 
ported here apply strictly for spherical accretion. The assump- 
tion of spherical accretion is not necessarily a bad one, since 
late-time accretion settles after many sound-crossing times, and 
the flow can recover from modest asymmetries arising in the ex- 
plosion. Furthermore, unlike accretion in a binary system, the 
accreting material does not acquire angular momentum from 
inspiral, and the specific angular momentum is limited to pre- 
explosion values, which may not be too large. On the other 
hand, late-time accretion is shaped in part by the luminosity at 
early times, for which disk-like rather than spherical flow may 
be more suitable (Mineshige et al. 1997). The radiative effi- 
ciency of hypercritical disk accretion is not well constrained by 
theory, and assessing its implications on the prospects of de- 
tecting black hole emergence is beyond the scope of this paper. 

3. MODEL RESULTS: LUMINOSITY AT EMERGENCE 

With the model of § 2 we proceed to calculate the time of 
emergence and corresponding luminosity as a function of pro- 
genitor mass, t BH (M„) and L acc {t B H) = L BH (M„), respec- 
tively. 

3.1. Emergence Luminosity with no Radioactive Isotopes 

The most favorable situation for observing black hole emer- 
gence, is when radioactive heating is negligible and emergence 
occurs once envelope recombination stratifies. The results for 
the accretion luminosity (Eq. [2]) based on the effective time 
tdust (Eq. [16]) and emergence time, t B H (Eq. [17]) as a func- 
tion of progenitor mass are shown in Figure 2. In the calcula- 
tions we varied the fraction of explosion energy spent on un- 
binding the star (/ = 0.25, 0.50 and 0.75), and the fraction of 
kinetic energy deposited in the helium layer (fn e = 0.01 and 
0.03). We set the progenitor radius at Rq = 5 x 10 13 cm and 
the fiducial fraction of thermal energy as f t h = 0.1. 

The accretion luminosities at emergence varies between sev- 
eral 10 35 to 3 x 10 37 ergs s _1 for the most optimistic models. 
The results are quite sensitive to both effective parameters / 
and /ne, and, for a given progenitor mass, can vary by more 
than an order of magnitude between extreme assumptions. In 
spite of these uncertainties, we can identify some global trends. 
First, for a given total explosion energy and black hole mass 
(given value of /), the luminosity at emergence increases with 
decreasing /h c : less energy to the helium layer means more 
bound mass which is also slowly moving as a reservoir for 
late-time accretion. Within the parameterization of the model, 

tdust oc / Ho 3/4 while t BH is very weakly dependent on / He , 

thereby yielding the relation L BB oc f^J 72 , which holds for 
progenitors of all masses. 

The luminosity at emergence also increases with increasing 
/, and that the sensitivity to this factor is at least as large as 
that to / Hc - Although a larger available energy for unbinding 
the star leads to a lower mass black hole, it implies a larger 
helium-rich layer in the ejecta and less kinetic energy for the 
ejecta (which scales as 1 — /). In the parameterization equa- 
tions (2) and (16-17) a thick and slow helium-rich layer (i.e., a 
larger reservoir of material for late-time accretion) outweighs 
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FIG. 2. — Parameter survey for the accretion luminosity at emergence, 
Lgu, as a function of progenitor mass, M r , for the model with no radioac- 
tive heating. Parameters surveyed are / and /h c , with / = 0.25, — dashed 
lines, / = 0.50 — thick solid lines, and / = 0.75 — thin solid lines, and with 
/ He = 0.01, and 0.03. 

the mass of the black hole in determining the magnitude of 
the accretion luminosity. For realistic values of / we find that 
roughly L B h oc /. This result reinforces previous expecta- 
tions (Zampieri et al. 1998a; Balberg et al. 2000) that super- 
novae with low inferred explosion energies (low expansion ve- 
locities) are the most promising in terms of observing black 
hole emergence. Note that since the actual values of black hole 
and helium-layer masses also depend on the progenitor mass, 
changing / does not just scale the L BB (M*) curve (as chang- 
ing /n e does) but also modifies the shape of this function. The 
dependence of emergence luminosity on progenitor mass can 
also be understood in terms of the competition between black 
hole mass and ejected helium layer mass. 

3.2. Emergence Luminosity with Radioactive Heating 

If radioactive decays dominate over accretion as the heating 
source after recombination, emergence can be significantly de- 
layed. We determine emergence time and luminosity by com- 
paring the evolution of the accretion luminosity (eqs. [2] and 
[16]) and the radioactive luminosity (eqs. [13] and [18]) for 
a given Ms6 Co (eq. [12]) in the ejecta. The emergence lumi- 
nosity as function of progenitor mass, L B h(M*), is shown 
for three combinations of / and /n e in Figure 3; these are 
{/ = 0.5 / He = 0.01},' {/ = 0.75 / He = 0.01} and 
{/ = 0.25 ,/ho = 0.03}, which are intermediate, most op- 
timistic and least optimistic, respectively, in terms of the ac- 
cretion luminosity they predict. In each case we compare four 
assumptions about the abundance of 56 Co : none (i.e., the re- 
sults in § 3.1), Ms6 Co = 2 x 10~ 3 M© for all progenitors, 
and a black hole mass-dependent abundance according to equa- 
tion (12), with 13 = 2 and 8 (labeled with the value of (3). 

The dominant feature of Figure 3 is that the qualitative effect 
of including radioactive heating is similar for all three cases 
(and, indeed, is similar for other combinations of / and /h 
not shown here). When a constant Ms6 Co ~ 2 x 10~ 3 M Q 
is assumed for all progenitor masses, the emergence luminos- 
ity is decreased by a factor of 5-50 with respect to its value in 
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FIG. 3. — Accretion luminosity at emergence, Lgu, as a function of progenitor mass M* for three combinations of / and /jj and different assumptions about 
the abundance of 56 Co in the ejected envelope. Curves are labeled with the mass of Ms6 Co if independent on M* and by the value of fi used for a M56 Co (M*) 
according to its value in equation (12). 



the radioactive-free case for the same {/, /h c } combination; 
a greater decrease corresponds to a lower value of Lbh in the 
radioactive-free case. This similarity results because the time 
of emergence is controlled by the exponential decline of the 
56 Co heating rate; the emergence times are then constrained to 
the rather narrow range of 1200 ± 200 days, hence constraining 
the value of L bh as well. 

When the amount of 56 Co is assumed to decline with in- 
creasing black hole mass, the Lbh curve ascends from the 
Mseco = 2 x 10~ 3 M Q values for lower mass progenitors 
to the radioactive-free values for higher mass progenitors; the 
slope of this ascent is, of course, determined in our parameter- 
ization by the value of [3. The net effect of such a black hole 
mass-dependent 56 Co abundance is to suppress the contribu- 
tion of the lower mass range of progenitors to observable black 
hole emergence events. The implications of this suppression 
on the rate of observable events of black hole emergence may 
be quite significant, if a realistic progenitor mass distribution is 
weighted towards the lower mass stars, these progenitors dom- 
inate the supernova rate. 

It is important to comment on the implications of 57 Co and 
44 Ti . If those isotopes are present with abundances which scale 
with M5 6Co similar to the ratios inferred in SN1987A (Woosley 
& Timmes 1996), they would dramatically alter the prospects 
for observing black hole emergence. It is unclear to what ex- 
tent such scaling is justified, since 57 Co and 44 Ti are synthe- 
sized even deeper than 56 Ni , so scaling is mostly useful for 
placing an upper limit on the 57 Co and 44 Ti abundances. If 
these abundances do scale according to M56 Co = 2 x 10 -3 M Q 
and the SN 1 987 A ratios, both become important for radioactive 
heating at about 1000 days after the explosion (Balberg et al. 
2000), with heating rates of L57 Co w La Ti < 10 35 ergs s _1 . 
Heating from 44 Ti decay (positron heating) should remain 
at this level for tens of years after the explosion. Thus, if 
Lbh estimated from comparison with Ls6 Co is slightly less 
than a several 10 35 ergs s _1 , unveiling the accretion lumi- 
nosity would be a more difficult task since the total luminos- 
ity will include additional components (rather than follow a 
simple power law decay). If Lbh is significantly less than 
10 35 ergs s _1 , as in the case of / = 0.25 / Hc = 0.03 and 
Mn6 Co = 2 x 10~ 3 M Q , the accretion luminosity would be 
sufficiently subdominant so that detection would be impossi- 
ble. Fortunately, we can avoid detailed investigation of the 
effects of 57 Co and 44 Ti since for HST capabilities an emer- 
gence luminosity below 10 35 ergs s _1 cannot be detected be- 



yond the Local Group (see below). Provided that scaling of 
57 Co and 44 Ti abundances as in SN1987A applies, any emer- 
gence luminosity that is low enough to be obscured by 57 Co and 
44 Ti heating is undetectable in any case. 

4. RATE OF OBSERVABLE BLACK HOLE EMERGENCE 

Having derived the function Lbh(M^), we now proceed to 
estimate the event rate. We define an "event" in this context as 
a case where the accretion luminosity at emergence is above the 
threshold of a given instrument (the total luminosity is twice as 
large, of course, but we assume conservatively that a factor of 
two margin would be needed to certify the nature of the accre- 
tion luminosity, especially if the power-law dependence is to be 
identified). The detector threshold translates into a minimum 
emergence luminosity which can be detected from any given 
distance. We will use the HST STIS camera capabilities as a 
our reference and require a threshold apparent bolometric mag- 
nitude of tub = 28.5. The function Lbh(M^) provides the 
progenitor masses which have emergence luminosity exceeding 
that minimum luminosity; their relative fraction out of all pro- 
genitors which make core collapse supernovae is determined by 
assuming a present day mass function (PDMF) of progenitors. 
In this section we estimate this fraction and compare it with the 
observed supernova rate in the local universe in completing an 
estimate of events of observable black hole emergence. 

4. 1 . Events per Supernova as a Function of Distance 

For simplicity, we assume a single power law PDMF for the 
progenitors of core collapse supernova, similar to the initial 
mass function (IMF) of the Milky Way. There is some observa- 
tional evidence the mass function of massive stars in the Milky 
Way does follow such a power law with no sign of dependence 
on metallicity (Massey & Thompson 1991). In this parameteri- 
zation, the fraction F(M) of progenitors with mass greater than 
M follows the relation 

— cx M-t ; F(M) = C CC M^ +1 ; (19) 

C cc = (Mp +1 MP+Y 1 , 

where the value of the constant C cc arises from stipulating that 
core-collapse supernovae occur in stars with progenitor masses 
between M x and M 2 . 

We also assume that all local galaxies which host core- 
collapse supernovae share a single power law 7. By changing 
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FIG. 4. — The fraction of Type II supernovae for which black hole emergence can be observed as a out of the all Type II supernovae as a function of distance for 
three combinations of / and / He and the various assumptions about the 56 Co abundance (see Fig. 3) and 7 = 2.7. 



7 over a reasonable range of values we can assess the likely un- 
certainties embedded in this assumption. The value of 7 is cru- 
cial in determining the black-hole-per-core-collapse-supernova 
fraction, which is essentially the fraction of progenitors which 
leave behind a remnant with M > 3 M Q . By following Fryer 

6 Kalogera (1999) who adopted the Scalo (1986) value of 

7 = 2.7 as a nominal choice and assumed core collapse super- 
nova occur for progenitors with masses in the range 8 — 40 M , 
we find that the fraction of supernovae which make black holes 
is about 0.15. 

The combination of the function Lbh(M^) and the PDMF 
of progenitors allows to estimate what is the fraction of black 
hole forming supernovae for which emergence is observable as 
a function of distance. It can be further expanded to estimate 
the fraction of emergence events out of all core-collapse super- 
novae (most of which produce neutron stars). In Figure 4 we 
show this fraction as a function of distance from the Milky Way, 
(for a threshold apparent magnitude of 28.5), assuming no dust 
extinction along the line of site. The three graphs correspond 
to the three combinations of / and /h c shown in Figure 3, with 
each case including the different assumptions concerning the 
abundance of 56 Co . The results are binned in 2 Mpc intervals. 

At very close distances, all supernovae which create black 
holes yield an observable emergence, so the "event" per super- 
nova rate is simply the fraction of supernovae which make black 
holes. The characteristic fraction of ~ 15% for 7 = 2.7 does 
have some dependence on the value of /, which determines the 
threshold progenitor mass for producing a black hole. As dis- 
tance is increased, for some progenitor masses the emergence 
luminosity falls below threshold and the event fraction (out of 
all supernovae) decreases. The significance of radioactive heat- 
ing is clearly evident in the details of each graph, where a larger 
M56C0 naturally delays the emergence of the black hole, and so 
the luminosity at emergence is lower and therefore is observ- 
able over smaller distances. The importance of the distinction 
between the two exponential cases, (3 = 2 and (3 = 8 is also 
evident. For (3 = 8 radioactive heating does not affect the max- 
imum values of Lbh so the dominant contribution to cases of 
observable emergence, especially at large distances, is similar 
for the radioactive-free case and the (3 — 8 case. On the other 
hand the (3 = 2 model results in reducing the maximum value 
of Lbh, thereby reducing the maximum distance, D max , of 
observable emergence. 

Most notable is the difference between typical distances over 
which emergence is observable for the different combinations 
of / and /ne- These vary from D max rts 14 Mpc for the 



worst case {/ = 0.25 fu e = 0.03} (and radioactive free) to 
D m ax ~ 54 Mpc for the best case {/ = 0.75 / He = 0.01}, 
which roughly amounts to a factor of about 60 in volume, and 
hence in the event rate. In general, we find that the combi- 
nation of the uncertainty in / and /ne and in the abundance 
of 56 Co places D max roughly in the range of 10 — 50 Mpc, 
corresponding to a span in the expected event rate of two or- 
ders of magnitude. In this last assessment we exclude the case 
of {/ = 0.25 ./hc = 0.03} and M56 Co = 2 x 10~ 3 M G 
for which black hole emergence is essentially unobservable be- 
yond the Local Group. Recall that this is also the combination 
where the accretion luminosity may be obscured if the amount 
of 57 Co and 44 Ti is not negligible. If this extreme combination 
is representative of black hole forming supernovae, the rate of 
observable emergence is very small. 

The slope of the PDMF, 7, has some influence over the frac- 
tion of black hole-forming supernovae and over the observable 
event rate. For values of 7 between 2 and 3, the fraction of 
core-collapse supernovae which form black holes varies be- 
tween 0.25 and 0.125, respectively. The point is that the value 
of 7 does not affect the maximum distance for observation for 
a given {/ /n e } model. Figure 5 compares the event per core 
collapse supernova fraction for {/ = 0.5 fn e = 0.01} for 
7 = 2, 2.7, and 3. In all cases the fraction at any given distance 
approximately follows the total fraction of progenitors which 
make black holes. Note that the graphs for 7 = 2 and 7 = 3 
are almost identical, after the y-axes are scaled 2: 1 .In total, the 
uncertainty in 7 accounts for a relatively modest uncertainty of 
a factor of 2 in the event rate. 

An additional modest source of uncertainty is the lower limit 
on progenitor masses which make core collapse supernovae. 
For 7 = 2.7, stars between 8 and 10 M Q make up about one 
third of all those between 8 — 40 M Q , so the fraction of su- 
pernovae from higher mass progenitors out of the total could 
be sensitive to our choice of the lowest mass progenitors by as 
much as 50%. On the other hand, unless the PDMF is excep- 
tionally flat, the exact value of the upper mass limit is not too 
important, since stars at the very high mass end are very rare. 

4.2. Eliminating Type Ib/c Supernovae 

So far we referred to estimates of the event rate per core- 
collapse supernova, which is the expected fate of massive pro- 
genitors. These include also Type Ib/c supernovae which are 
not good candidates for observing black hole emergence (Bal- 
berg et al. 2000): if the hydrogen envelope has been ejected 
prior to explosion, practically all the kinetic energy is deposited 
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FIG. 5. — The fraction of Type II supernovae for which black hole emergence can be observed as a out of the all Type II supernovae as a function of distance for 
for 7 = 2.7 (nominal) and 2 and 3, the various assumptions about the 56 Co abundance (see Fig. 3), and {/ = 0.5 /h c = 0.01}. 



in the helium layer (/ He = 1)- Effectively, the evolutionary 
option of a Type Ib/c supernova reduces the fraction of core 
collapse supernovae in which black hole emergence could be 
observable. 

Studies of single and binary stellar evolution indicate that 
the likelihood that sufficient mass loss will occur during the life 
time of a massive star depends mainly on the presence and prop- 
erties of a binary companion; the mass of the star appears to 
play a secondary role. In our analysis we use the simplest clas- 
sification offered by Willstein & Langer (1999), where massive 
stars in binaries lose all their hydrogen envelope and result in a 
Type Ib/Ic supernova, and single massive stars do not lose all of 
the envelope, giving rise to a SNII. This simplest interpretation 
is that the probability of a massive star yielding a SNII, P/j, is 
independent of its mass. With this assumption, the calculations 
of the fractions shown in Figures 4-5 are then the fractions per 
SNII, related to the fraction per core-collapse supernovae by a 
factor of Pji. In the following we will set Pjj = 0.85, which is 
roughly consistent with theoretical estimates of supernova rates 
(Thielemann 1994; Cappellaro et al. 1997). 

4.3. Rate of Events in the Local Universe 

If the efficiency of discovering supernovae in the local uni- 
verse (D w 50 Mpc) were perfect, event-per-supernova frac- 
tions such as those shown in Figure 4 should be compared to 
the estimated supernova rate per galaxy. These are somewhat 
constrained based on observation and theory, and generally ex- 
pressed in units of SNu (ISNu = supernovae per century per 
10 10 Lq). The estimated rates (Thielemann 1994; Cappellaro 
et al. 1997) are dependent upon galaxy type, where core col- 
lapse supernovae appear to be limited to spiral and irregular 
galaxies. The combined theoretical rates (in SNu) for spiral 
galaxies are: SNell - about 0.7 ± 0.3, and Type Ib/c - about 
0.14 ± 0.07 (with some dependence on specific galaxy type). 
In irregular and peculiar galaxies there appears to be a slightly 
larger fraction of Type Ib/c supernovae (SNelb/Ic), whereas in 
elliptical there are hardly any core collapse supernovae at all. 

However, the actual constraint on observing black hole emer- 
gence rate is the rate at which supernovae are discovered, rather 
than the rate at which they occur. Clearly, any attempt to un- 
cover the emergence of the black hole would be conditional 
on the supernova being found in the first place. The efficiency 
of discovering supernovae is imperfect due to both incomplete 
monitoring of the sky and line of sight limitations (background, 
extinction etc.). For example, unfavorable inclination of the 
host galaxy appears to cause a deficiency of observed super- 



nova in galactic cores with respect to those found which are 
offset from the galactic centers (Cappellaro et al. 1997). 

4.3.1. Rate of Observed Supernovae 

A recent compilation of supernovae catalogs (Barbon et 
al. 1999) suggests that the recent (1997-1998) rate of ob- 
served supernovae identified as core collapse in the local uni- 
verse is roughly 10~ 4 Mpc~ 3 yr _1 although the uncertainty in 
this estimate is quite large. This number includes both SNell 
and SNelb/Ic, with a ratio of about 5:1, respectively. Within 
the statistical uncertainties of the search, these rates appears 
to be independent of distance up to 50Mpc. Such distance- 
independence arises from the limits of supernovae searches (in- 
complete coverage of the sky, methods of the searches), which 
currently dominate over the dimming of the more distant su- 
pernovae. A rate of 10~ 4 Mpc~ 3 yr _1 is roughly 10-20% of 
the total rate of core collapse supernovae which are expected to 
occur in the local universe. 

While we will use this estimate of 10~ 4 Mpc~ 3 yr _1 core 
collapse supernova as a guideline in our calculation, it is clear 
that at very small distances any approximation of a uniform ob- 
served supernova rate must break down. In particular, the ac- 
tual supernova rate up to several Mpc from the Milky Way is 
dominated by a few individual galaxies: most notably, the star 
burst galaxies NGC 253 (3.0 Mpc) and M82 (3.5Mpc). We rely 
on the estimates of Becklin (1991), who evaluated the core- 
collapse supernova rate up to 12 Mpc from the properties of 
observed galaxies. His estimate was based on a rather high 
theoretical rate of core collapse supernovae: 6 per century per 
10 10 Lq in IR (the contribution of massive stars to IR emission 
exceeds that in optical). Recent surveys seem to favor a rate 
which is about one half of that. We therefore modify those esti- 
mates by a factor of 0.5. The rates of core collapse supernovae 
that we use for this range are: 2-4 Mpc - 0.3yr _1 ; 4-6 Mpc - 
0.02yr- 1 ; 6-8 Mpc - O^yr" 1 ; 8-10 Mpc O.lSyr" 1 ; and 10-12 
Mpc - 0.3yr _1 . This last value fits onto the estimated observed 
rate of 10~ 4 Mpc~ 3 yr _1 assumed to hold beyond these clos- 
est distances. We do not include here supernovae in the Local 
Group (0-2 Mpc), since they are dominated by the Milky Way 
and for the most part would be unobservable. 

4.3.2. The Rate of Observable Black Hole Emergence Events 

We first assume naively that the probability of detection of 
a SNII is independent of progenitor mass (hereafter the UNBI- 
ASED model). The event fraction per occurring core collapse 
supernova in § 4. 1 then translates simply into an event fraction 
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per observed core collapse supernovae. Using the estimates 
for the observed rate above, we then calculate the total rate 
of observable black hole emergence events, including a factor 
of 0.85 to count only SNell where emergence is observable. 
These estimated rates as a function of /, /h and the assump- 
tions regarding the 56 Co abundance are shown in Table 1. 

Table 1 

Predicted rates (yr _1 ) for observable black hole 
emergence in the "naive" model (see text). 
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The results in Table 1 indicate that the likely rate of detect- 
ing of black hole emergence in supernovae in nearby galaxies 
clearly depends on all three of our main unknowns - the fraction 
of energy spent on unbinding the stars, distribution of kinetic 
energy between the hydrogen envelope and the helium layer, 
and the abundance of 56 Co . Within the limits we used in our 
models, each parameter alone imposes an uncertainty of a fac- 
tor of a few in the event rate, whereas the combined effect of all 
three accumulates to about two orders of magnitude of uncer- 
tainty. The prediction for the event rate varies between as many 
as a few events per year to a few events per century. In the for- 
mer case, one expects that several of the SNell discovered each 
year will lead to an observable black hole emergence, were they 
to be monitored for months after the explosion. Note that for 
most of the parameter range, the event rate should at least be 
one event every few years; rates lower than that are confined al- 
most uniquely to the worst case model for radioactive heating, 
when a constant Mse Co = 2x 10~ 3 M Q for all supernovae is 
assumed. 

4.3.3. Effects of Bias in Observability of Supernovae and 
Emergence 

The estimates above may be too naive in their prediction of 
the rate of observable events of black hole emergence, since it 
explicitly assumes that all SNell are equally likely to be de- 
tected. A more realistic point of view is that the probability of 
detecting SNell is not uniform, but rather is dependent on its 
peak luminosity. Recall that the best known candidate for ob- 
serving black hole emergence, SN1997D, is considered signifi- 
cantly subluminous in comparison with typical SNell. While it 
is unlikely that underluminous supernova dominate the total su- 
pernova rate (Cappellaro et al. 1997 limit them at 20 - 30% of 
the total rate), they may represent a large fraction of the black 
hole forming explosions. If, on average, black hole forming 
supernovae have a lower probability of being discovered, the 
fraction of events per observed SNII is obviously reduced with 
respect to our first estimate. 

The luminosity of supernovae during the diffusion and re- 
combination phase is known to depend on the explosion energy 
and progenitor radius and structure (Arnett 1996), and indeed, 
observations indicate a heterogeneous luminosity distribution. 
In particular, relatively compact progenitors will result in an 
underluminous light curve even if the explosion energy was 
quite large, as was the case in the very well studied SN1987A. 
Nonetheless, it is reasonable to expect some correlation in the 



detectability of supernovae with explosion energy. For lack 
of a better model, we assume that the probability of detecting 
the supernova is proportional its typical luminosity. We further 
assume that this typical luminosity is proportional to explosion 
energy deposited in the envelope (i.e., to (1 — f)E exp in our 
parameterization). This approximation is roughly correct for a 
given progenitor (Arnett 1996; Balberg et al. 2000). We show 
the quantitative effect of such a bias against detecting low lu- 
minosity supernovae on the resulting fractions of event per ob- 
served SNII as a function of distance in Figure 6b, for 7 = 2.7, 
{/ = 0.50 /ho = 0.01} and after normalizing the fraction of 
all observed supernovae in the range 8 — 40 M Q to unity. The 
immediate result in this model (hereafter the E-BIAS model) is 
a significant reduction of the typical fraction of event per ob- 
served supernova, which is now down to about 0.02 instead of 
0.15. Also, the relative impact of the bias is larger in the cases 
with a progenitor mass dependence of the 56 Co abundance. 
This is quite expected, since in these cases the effect of 56 Co is 
limited to the lower mass progenitors, so the relative contri- 
bution of the higher mass progenitors in the event is greater. 
The corresponding predicted rates under this E-BIAS model 
for the various {/ / H o} combinations are shown in Table 2. 

Table 2 

predicted rates (yr _1 ) for observable black hole 
emergence in the model with an explosion energy bias 
(see text - model e-bias). 
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Overall, we find that biasing the supernova detection prob- 
ability according to explosion energy has a profound effect on 
the expected rate of observable black hole emergence. In gen- 
eral, we find that the event rate drops by a factor of 5-15 with 
respect to an unbiased detection probability, as depicted in Ta- 
ble 1 . Note that the effect is more pronounced for larger values 
of /, because the luminosity is driven by the energy that is not 
spent on unbinding the star, i.e., on 1 — /. The typical event 
rates with these revised assumption are mostly once per several 
years for /n e = 0.01 and once per tens of years for /h c = 0.03. 
Once again, a constant Ms6 Co = 2 x 10~ 3 M Q for all progen- 
itors would be especially destructive in terms of decreasing the 
potential event rate. 

An additional source of potential bias that would reduce the 
observability of black hole emergence is luminosity due to in- 
teraction of the ejected envelope with circumstellar material 
(CSM). If the CSM is dense enough, a significant luminosity, 
Lcsm, will arise as the shock wave propagates through it. In- 
teraction with the CSM is the preferred model for the subclass 
SNelln, which are characterized by a high, slowly decaying, lu- 
minosity after the recombination peak. In particular, Chevalier 
& Fransson (1994) demonstrated that a dense CSM could sus- 
tain a luminosity at a level of Lqsm ~ 10 40 ergs s _1 for sev- 
eral hundred years after the explosion. A luminosity at this level 
would black hole would obscure any potential accretion lumi- 
nosity. Indeed, this was quite possibly the case for SN1994W 
(Sollerman et al. 1998), which otherwise could have been a 
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FIG. 6. — The fraction of observable Type II supernovae for which black hole emergence can be observed as a function of distance out the total number of 
observable Type II supernovae for three sets of assumptions about supernova observability and accretion luminosity. These are (from left to right) - UNBIASED: 
(same as left panel in Fig. 4); E-BIAS: underdetection of low energy explosions, and E+M-BIAS: underdetection of low energy explosions combined with CSM 
obstruction of black hole emergence (see text). In all cases / = 0.5, /h c = 0.01 and 7 = 2.7; all variations for the 56 Co abundance are shown (see Fig. 3). 



promising candidate for observing black hole emergence. Cap- 
pellaro et al. (1997) estimated that Type Iln supernovae are 
relatively rare, and probably make up only 2 — 5% of the total 
SNell rate. However, they note that owing to the relatively high 
luminosities of these supernovae, in recent years their contribu- 
tion to the observed rate of SNell is 15 — 20%; their influence 
over the present study must therefore be examined. Moreover, 
since the CSM is essentially the residue of mass loss from the 
progenitor during the red giant phase, more massive progenitors 
are more likely to have produced denser CSM. 

While the rate of mass loss in the star's wind, M W i n d, is ex- 
pected to increase with progenitor mass, many other factors will 
contribute to the overall properties of the wind and the resulting 
CSM. We note that there was no evidence of CSM interaction at 
any stage in the observation of SN1997D (Benetti et al. 2001). 
For the purpose of a quantitative estimate of the potential in- 
fluence of CSM interaction on the prospects of observing black 
hole emergence, we assume that the probability of the accretion 
luminosity not being obscured by a high Lqsm decreases for 
larger progenitors and scales as M~ e . We calibrate 9 by im- 
posing the following combination: assume that (a) significant 
Lcsm occurs in 5% of SNell and that (b) if all supernovae with 
CSM interaction are observed, they make up 20% of observed 
SNell. Imposing these assumptions on the E-BIAS model de- 
scribed above, we find that for 7 = 2.7 this combination is best 
reproduced with « 0.5. By including the effect of possible 
CSM interaction on the light curve according to this recipe and 
renormalizing the overall observed SNell rate, we arrive at yet 
another set of event per observed SNell fractions as a function 
of distance (hereafter the E+M-BIAS model). These sets are 
presented in Figure 6c. Note that these new fractions are natu- 
rally lower than the corresponding values in the E-BIAS model. 
The total predicted rates under these assumptions are given Ta- 
ble 3. 

Quantitatively, the potential of CSM interactions contributes 
an additional effect of reducing the predicted event rates by 
~ 30% with respect to Table 2. It also further reduces the event 
rate in the (3 = 8 model with respect to the radioactive free 
one, emphasizing the significance of even a small abundance 
of 56 Co in the lower mass progenitors. While the additional 
mass dependent bias is not as important as the explosion en- 
ergy bias, owing to the different functional dependence on M„, 
the effect of CSM interactions - as parameterized in our model - 
reduces the rate of observable black hole emergence to, at best, 
one per four years. For most of a parameter range, we conclude 



that with using HST capabilities and current limits of supernova 
searches this rate is probably not better than once per decade. 

Table 3 

Predicted rates (yr- 1 ) for observable black hole 
emergence in the model with an explosion energy and 
wind bias (see text - model e+m-bias). 
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5. CONCLUSIONS AND DISCUSSION 

In this work we attempt a first estimate of the rate at which 
black hole emergence in supernovae may be observed. We 
used an analytic model and a parameter survey to examine the 
competition between the expected accretion luminosity and the 
other sources of power for the supernova light curve, most no- 
tably radioactive heating. In general, we show that the event 
rate is larger if the relative fraction of the total kinetic en- 
ergy driving the expansion (not spent on unbinding the star) 
is smaller and if the fraction of this kinetic energy deposited in 
the helium layer (which is the source of late-time accretion) is 
also smaller. Furthermore, the event rate is larger for reduced 
amounts of radioactive heating in the envelope. 

With the efficiency of current supernova searches and the ca- 
pabilities of the HST for detecting the accretion luminosity at 
emergence, we find that under optimistic assumptions, the rate 
of "events" - successful observations of black hole emergence 
- may be as large as a few events per year. If Type II super- 
novae in current searches are detected with no dependence on 
progenitor mass and explosion energy, then for most of our pa- 
rameter range, the event rate should be better than once per 
two years. However, in the more likely event that detection 
of supernovae in current searches is weighted towards brighter 
supernovae, this rate is dramatically reduced. The more mas- 
sive progenitors which lead to black hole formation are also 
likely to produce less luminous explosions, and would there- 
fore be underrepresented in present searches. The best known 
potential candidate for allowing observable black hole emer- 
gence, SN1997D, was such a subluminous explosion. If a bias 
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against discovering supernovae from explosions of more mas- 
sive progenitors exists, it reduces the potential rate of finding 
black hole emergence in supernovae as well. We estimate that 
the event rate in the presence of such a bias is not larger than 
once per few years. We also find that some additional bias may 
arise due to obstruction of the accretion luminosity by inter- 
action of the expanding envelope with circumstellar material. 
Such obstruction is more likely for the more massive progeni- 
tors which experience greater winds during their evolution; we 
view SN1994W as a likely candidate for potential emergence 
that was obscured by luminosity due to circumstellar interac- 
tion. We find that if probability of significant contribution to 
the light curve from interaction of the ejecta with circumstellar 
material is larger for more massive progenitors, the fraction of 
observable black hole emergence per observed Type II super- 
novae is reduced by an additional factor of 0.3 — 0.7. In this 
later case the actual event rate for present detection capabilities 
and limitations of supernovae searches is unlikely to be greater 
than one per decade. 

Our findings suggest that, at present, the prospects of ob- 
serving black hole emergence in Type II supernovae are small, 
but not negligible. The immediate consequence is that if a low 
energy Type II supernova is discovered, high priority must be 
given to following it for several months until beyond recom- 
bination. If a small (~ 10~ 3 M Q or less) amount of 56 Co is 
then identified in the tail of the light curve, it would be prudent 
to schedule HST observations of the declining light curve with 
hope of observing the emergence of the black hole. 

These rates might be somewhat increased by detailed spec- 
tral analyses, although a successful spectral analysis requires 
a larger-than-threshold total luminosity. Since luminosity due 
to radioactive heating and CSM interaction is generated in the 
outer (and cooler) part of the expanding envelope, their spectra 
tend to be mostly in the R-band. Furthermore, in typical Type 
II supernovae the light curve is observed to become redder at 
later times, as the effective temperature of the emitting region 
decreases. By contrast, the accretion luminosity is generated 
deep in the ejecta, so its effective temperature should be that 
of partially ionized envelope material - 5000-10000 degrees. 
Frequency-dependent calculations of the accretion luminosity 
have not been carried out, but the intensity should peak in the 
visible or near IR. A light curve in which the R-band luminosity 
follows a decline rate expected from 56 Co decay while the vis- 
ible band luminosity appears to be decline more slowly could 



indicate the presence of the black hole even if accretion is still 
subdominant in the total bolometric luminosity. 

Perhaps the most important conclusion of our analysis is that 
detectability of the supernovae themselves appears to be the 
dominant factor in the rate of observing black hole emergence, 
if indeed black hole producing supernovae are fainter than av- 
erage. Clearly, a search for black hole emergence will benefit 
greatly from any dedicated faint supernova project, indirectly 
also from one aimed at high red shift explosions. If the bias 
we assumed against discovery of fainter supernovae (Table 2) 
can be at least partially lifted by such a search, we may be able 
to approach the higher rates found for an unbiased survey (Ta- 
ble 1). An additional boost to the event rate would arise with a 
next generation instrument with improved capabilities. While 
57 Co and 44 Ti abundances may not allow detection of emer- 
gence luminosities lower than ~ 10 35 ergs s _1 , a more sen- 
sitive instrument will certainly increase the range over which 
higher emergence luminosities can be observed. In particu- 
lar, we find that if the threshold apparent magnitude of the in- 
strument can be improved to m v ~ 30 (possible threshold for 
NGST), the expected event rate is increased by a factor of 6-8. 

In conclusion, we suggest that observing emergence of a 
black hole in Type II supernovae is potentially possible even 
with present capabilities. In view of the relatively low rate, high 
priority should be assigned if a likely candidate supernova (low 
energy, low 56 Co ) is discovered. The event rate would sig- 
nificantly enhanced with future improved instruments, and also 
through dedicated nearby faint supernova searches with present 
capabilities. Such a search based on ground instruments is cur- 
rently being initiated (M. Turrato 2000, private communica- 
tion), and could also be included as part of other projects aimed 
at observing high redshift supernova. In view of the importance 
of achieving direct confirmation of the supernova - black hole 
connection, we believe such efforts should be encouraged. 
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